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In the last few years, radio detection of cosmic ray air showers has experienced a true renaissance, becoming 
manifest in a number of new experiments and simulation efforts. In particular, the LOPES project has success- 
fully implemented modern interferometric methods to measure the radio emission from extensive air showers. 
LOPES has confirmed that the emission is coherent and of geomagnetic origin, as expected by the geosynchrotron 
mechanism, and has demonstrated that a large scale application of the radio technique has great potential to 
complement current measurements of ultra-high energy cosmic rays. We describe the current status, most recent 
results and open questions regarding radio detection of cosmic rays and give an overview of ongoing research and 
development for an application of the radio technique in the framework of the Pierre Auger Observatory. 



1. Introduction 

About 40 years ago, Jelley et al. measured 
pulsed radio emission originating from extensive 
air showers (EAS) for the first time. As the ra- 
dio technique proved to be too difficult to handle 
with the technical limitations of the 1960s and 
1970s, however, the interest in radio detection of 
cosmic rays diminished completely within the fol- 
lowing decade. It was only recently, now having 
powerful digital technology at our disposal, that 
the concept of radio detection of cosmic rays ex- 
perienced its renaissance |2). By now a number 
of new projects dedicated to the measurement 
of radio emission from EAS, most prominently 
the LOPES project 02] and the CODALEMA 
project [5], have been established. 

The radio technique for measuring cosmic rays 
has a number of merits in its own right: it can 
watch for EAS with nearly 100% duty cycle, even 
in populated areas, and it measures a bolomet- 
ric signal that is only very slightly attenuated in 
the atmosphere, thus allowing the observation of 
highly inclined showers. But naturally, the most 
interesting application is to combine the tech- 
nique with other detection methods, in particu- 
lar ground-based particle detector arrays and air 
fluorescence measurements. Each of these tech- 
niques yields different observables, and a combi- 
nation of the methods allows so-called "hybrid 
detection" of cosmic rays, yielding much more in- 
formation than the individual techniques alone. 

In this article, we review the goals, the status 
and the results so far gathered within the LOPES 
project followed by an outlook on the application 
of the radio technique on large scales for the mea- 
surement of ultra-high energy cosmic rays. 

2. The LOPES project 

The LOPES project was initiated in 2001, em- 
ploying a LOFAR 1 Prototype Station to detect 
radio emission from EAS in coincidence with the 
well-established KASCADE particle detector 
experiment. The goals of the LOPES project are 
(i) to deliver the "proof of principle" for radio de- 

1 LOFAR is a revolutionary digital radio interferometer for 
the 10-200 MHz range being set up in the Netherlands. 



tection of cosmic rays with modern interferomet- 
ric methods, (ii) to study and calibrate the radio 
emission in the energy regime up to ~ 10 18 eV, 
and (iii) to develop and optimise the radio tech- 
nique for large scale application at ultra-high en- 
ergies. To achieve these goals the LOPES project 
combines both experimental activities and dedi- 
cated theoretical efforts. 

2.1. Simulations 

Historical studies of the theory of radio emis- 
sion from EAS already pointed to a geomagnetic 
process as the dominant emission mechanism . 
These studies, however, were not detailed enough 
to serve as a basis for the interpretation of LOPES 
experimental data. We therefore started our own 
simulations of radio emission from EAS, concen- 
trating on the dominating geomagnetic emission. 
The scheme in which this mechanism is described 
is that of "coherent geosynchrotron radiation" [2] 
from shower electrons and positrons that are de- 
flected in the earth's magnetic field. 

In a first step, we studied the emission with a 
frequency-domain analytical model, analysing in 
particular the important coherence effects arising 
because the length scales present in the air shower 
are of the same order as the observing wavelength 
in the frequency range of LOPES |Sj. In a sec- 
ond step, we developed a detailed time-domain 
Monte Carlo simulation of the emission process [H] 
predicting the radio emission and its dependence 
on different air shower parameters with unprece- 
dented detail [TU] . 

Important predictions of these simulations are 
that (i) the radio signal should scale approx- 
imately linearly with the primary particle en- 
ergy in the coherent frequency range (see Fig. 
0, (ii) the frequency spectrum declines rather 
steeply to high frequencies, making low frequen- 
cies favourable for detection (see Fig.|2J), (iii) the 
electric field strength should decrease exponen- 
tially with radial distance to the shower centre, 

(iv) the dependence on the strength and geometry 
of the magnetic field is much more subtle than in- 
tuitively expected and mostly shows up in the po- 
larisation characteristics of the radio signal, and 

(v) inclined showers exhibit a much larger radio 
"footprint" , making them especially favourable 
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Figure 1. Scaling of the field strength emitted at 
2/rj = 10 MHz by a vertical EAS as a function of 
primary particle energy |1U|. From top to bottom: 
20 m, 100 m, 180 m, 300 m and 500 m to the 
north from the shower centre. The field strength 
scales as a power-law with index close to unity, 
as expected for coherent emission. 

for radio detection (see Fig. OJ). The essence of 
these elaborate simulations is summarised by a 
parametrisation formula that can be a useful tool 
to quickly estimate the strength of radio emission 
for a particular situation jl()j . 

Recently, we have developed a next-generation 
Monte Carlo code, replacing the parametrised 
description of air shower characteristics with a 
highly realistic CORSIKA based air shower 
model. First results are presented in |12j . 

2.2. Experiment 

LOPES uses so-called "inverted V" dipole an- 
tennas, each on their own sensitive to a large frac- 
tion of the sky, to measure radio signals in the 40- 
80 MHz band. A readout of the LOPES antennas 
is triggered whenever the KASCADE array regis- 
ters a large (approx. > 10 16 eV) air shower. This 
coincident measurement of the particle and radio 
components of an EAS is imperative for the un- 
ambiguous association of radio detections to air 
showers and, consequently, air shower parameters 
on an event-to-event basis. 

Using the air shower geometry as reconstructed 
by KASCADE, the radio data gathered by the 
individual LOPES antennas can then be timc- 
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Figure 2. Frequency spectra of a 10 17 eV ver- 
tical air shower at various distances to the north 
from the shower centre ^TJ|: From top to bottom: 
20 m, 140 m, 260 m, 380 m and 500 m. "Noise" in 
the incoherent high-frequency regime stems from 
the simplified air shower model and statistics. 

shifted and correlated appropriately in an offline 
analysis to let the array of antennas work as an 
interferometer looking into the specific direction 
of the air shower. This process is called "beam- 
forming" and is essential for increasing the signal- 
to-noise ratio of the radio measurements. For de- 
tails of the technical aspects of LOPES and data 
analysis procedures we refer the reader to 2J. 

In the initial LOPES setup, 10 linearly east- 
west polarised antennas were situated inside the 
KASCADE array of the KASCADE-Grande ex- 
periment ^21 ( re d triangles in Fig.QJ. This con- 
figuration, hereafter called LOPES10, collected 7 
effective months of data in the year 2004. The 
results derived from these measurements are pre- 
sented in section 1231 

After the measurements with 10 antennas had 
gained good statistics, the array was rearranged 
and extended to 30 linearly east-west polarised 
antennas situated mostly inside the KASCADE 
array (black triangles in Fig. [3J|. This configu- 
ration, hereafter called LOPES30, offers larger 
baselines (e.g., for a better angular resolution), a 
larger collecting area, and allows analyses with in- 
dependent "sub-arrays" of antennas, e.g., to mea- 
sure the lateral profile of the radio emission on a 
per-event basis. The most important enhance- 
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Figure 3. 10 MHz radio footprints of a verti- 
cal 10 17 eV air shower (left) and a 45° inclined 
10 17 eV air shower coming from the south (right) 
The signal increases by 0.25 fxV m _1 MHz -1 
per contour line. Inclined showers have a much 
larger radio footprint, making them particularly 
suitable for radio detection. 

ment from LOPES10 to LOPES30 was, however, 
that the LOPES30 array has been fully calibrated 
with an external reference source, allowing to fi- 
nally address the 40-year-old question of the ab- 
solute field strength of radio emission from EAS, 
in comparison with the simulation predictions. In 
addition, a detailed environmental monitoring al- 
lows to correlate the radio data to many external 
variables, including the static atmospheric elec- 
tric field at ground level [1*4] . 

Currently, LOPES is being reconfigured to 
make dual-polarisation measurements. These will 
allow a much better analysis of the angular corre- 
lations in the radio data, and a clear verification 
of the geomagnetic emission process. An addi- 
tional triggering by the Grande stations of the 
KASCADE-Grande experiment is also planned. 

2.3. LOPES10 Results 

The first analysis of LOPES10 data concen- 
trated on a limited number of very energetic 
events, selected by cuts on high (truncated) muon 
number and high electron number as provided by 
KASCADE. The results of this analysis ["""" were 
that (i) there is an unambiguous association of 
pulsed radio emission with air showers, and the 
KASCADE-reconstructed direction is consistent 
with the source of the radio emission, (ii) the ra- 
dio field strength correlates very well with the 
muon number, which in turn correlates well with 
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Figure 4. Layout of the LOPES10 and LOPES30 
array configurations relative to the KASCADE- 
Grande surface detector components. 

the primary particle energy, and (iii) the radio 
field strength correlates clearly with the angle 
between shower axis and earth's magnetic field 
(hereafter called geomagnetic angle). With this 
analysis, the "proof of principle" for radio detec- 
tion of EAS with modern interferometric meth- 
ods had been achieved. Additionally, the ob- 
served angular correlation with the geomagnetic 
field proved that a large fraction of the emission 
must be produced by a geomagnetic effect. 

A repetition of this analysis with much higher 
statistics and for multiple data sets selected with 
different cuts was repeated in ""£]. This elabo- 
rate analysis confirmed the earlier findings and 
constrained the correlations, e.g. on geomagnetic 
angle (see Fig. [5J), much more precisely. The re- 
sults have been summarised in a parametrisation 
formula describing the measured radio emission's 
dependence on the major air shower parameters. 
The scaling with primary particle energy is con- 
firmed to be almost linear (see Fig.[SJ|, as expected 
from theory for coherent emission. Also, the elec- 
tric field turns out to decrease exponentially with 
radial distance from the shower axis, again as ex- 
pected from the simulations. A direct compar- 
ison of the angular correlations with the simu- 
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Figure 5. Correlation of the radio pulse field 
strength (normalised by the assumed correlation 
on truncated muon number) with the angle be- 
tween shower axis and geomagnetic field 4 . 
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Figure 6. Correlation of the radio pulse field 
strength (normalised by the assumed correlation 
on geomagnetic angle) with the primary particle 
energy as reconstructed by KASCADE [3]. 

lations is difficult, however, because LOPESfO 
only measured the east-west linear polarisation 
component, thus missing a significant fraction of 
the emission depending on the air shower geome- 
try. Also, the field strength of LOPES10 data is 
not absolute calibrated, a problem which will be 
solved with LOPES30 data. 

Both experimental data and simulations show 
similarities, but also differences, to the parametri- 
sation of the historical experimental data given 
in 0. As the historical measurements were done 
with a very different setup and the documenta- 
tion of the data is not in all cases unambiguous, 
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Figure 7. Radio detection efficiency as a func- 
tion of mean radial distance between shower core 
and antennas as well as primary particle energy 
(left) and primary particle energy alone (right) 
[T^] . Detection efficiencies are for purely east- 
west polarised antennas. 

the differences are, however, not surprising. 

Another analysis has been carried out using the 
Grande stations (instead of the KASCADE ar- 
ray) to reconstruct the air shower parameters |16| . 
The air showers selected for this analysis have 
their core on average at much larger distances to 
the radio antennas than the ones selected in the 
previously described analyses. The most impor- 
tant results of this analysis are (i) LOPES does 
detect radio emission from EAS at distances as 
large as 500 m from the shower core for energies 
well below 10 18 eV, an important result for a pos- 
sible use of the radio technique on large scales, 
(ii) the radio pulse field strength depends signif- 
icantly on the accuracy of the input parameters 
(core position, shower axis), i.e., radio measure- 
ments have the potential to constrain these pa- 
rameters with high precision and thus help in the 
event reconstruction of classical particle detector 
arrays, (iii) an exponential decrease of the radio 
signal with radial distance to the core seems to 
describe the data also to large distances, and (iv) 
LOPES reaches a detection efficiency of approx. 
> 50% starting at approx. fO 17 eV (see Fig. 0, 
which is encouraging when taking into account 
that LOPESfO misses up to half of the emission 
because it measures only the east-west linear po- 
larisation of the radiation. 

As predicted by the simulations, inclined air 
showers pose a particularly interesting target for 
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Figure 8. Inclined air showers distribution in 
zenith angle and geomagnetic angle [TJ\. Blue 
crosses: events fully reconstructed by KAS- 
CADE, green triangles: events only geometrically 
reconstructed by KASCADE, red rhombs: events 
detected in radio. No showers from 90° and 270° 
azimuth angles (denoting east and west) are de- 
tected in radio. 

radio observations. On the one hand, their radio 
footprint is expected to be large, and on the other 
hand, inclined showers provide a larger "lever 
arm" to analyse the radio field strength correla- 
tions on zenith, azimuth and geomagnetic angles. 
Also, nearly horizontal air showers can be induced 
by neutrino interactions in the atmosphere, mak- 
ing them a particularly interesting field of re- 
search. In ^7] we have analysed a selection of 
air showers with zenith angles over 50°, an an- 
gular range where flat particle detectors such as 
those of KASCADE become inefficient. The anal- 
ysis has shown that (i) LOPES does detect very 
inclined air showers up to approx. 80° zenith an- 
gle, (ii) the radio detection efficiency rises with 
increasing zenith angle, and (iii) angular effects 
are indeed more pronounced in this zenith angle 
range. Figure |H1 shows the KASCADE triggered 
and radio detected events in the plane of azimuth 
angle versus geomagnetic angle. While statistics 
are still somewhat low, one can see two effects: 
First, there is a pronounced north-south asym- 
metry (more radio detections for showers com- 
ing from the south with large geomagnetic angle 
than from the north with small geomagnetic an- 
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Figure 9. Enhancement of the radio emission 
during thunderstorm conditions in comparison 
with events with similar air shower parameters 
recorded during fair weather |18j . 

gle). Such an asymmetry was also observed by 
the historical experiments and can be explained 
by the simulations when threshold effects play 
a role, as is certainly the case for LOPES. Sec- 
ond, there is a pronounced gap with no radio 
detections at azimuth angles of approx. 90° and 
270°, i.e., there are no radio detections of show- 
ers coming from east or west. As LOPES10 only 
measured the east-west linear polarisation of the 
emission, this is consistent with the simulation 
prediction that showers coming from east or west 
should be mostly north-south polarised. Addi- 
tionally, antenna gain effects could play a role. 

An analysis of radio events measured during 
thunderstorm conditions j!8| revealed that ra- 
dio emission from EAS seems to be amplified by 
the strong (of order kV/cm) atmospheric elec- 
tric fields present in thunderstorm clouds. Fig- 
ure El compares the field strengths of events dur- 
ing thunderstorms with that of events with sim- 
ilar air shower parameters measured during fair 
weather ("twin events"). In a number of cases, 
a clear enhancement of the radio emission is vis- 
ible. It also becomes clear that amplification of 
radio emission only occurs when very high elec- 
tric field strengths are present in the atmosphere. 
The fair weather electric field (of order 10 V/m) 
does not influence the radio emission from EAS. 
These effects are currently investigated in more 
detail using extensive simulations and LOPES30 
data. 
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3. The current understanding 

As of today, some properties of radio emission 
from EAS are well understood: The emission is 
clearly coherent as visible from the approximately 
linear scaling of field strength with energy seen 
in both experimental data and simulations |4I15| . 
Also, the fact that there is a pronounced corre- 
lation between the radio field strength and the 
geomagnetic angle indicates that a geomagnetic 
emission mechanism must be responsible for the 
dominant fraction of the emission |4I15I17| . 

Some theoretical predictions seem to be con- 
firmed but need to be investigated further with 
additional data: The radial dependence of the 
radio field strength is described well by an expo- 
nential decrease both in experimental data and 
simulations (7I10I16ITT7) . Also, as expected from 
theory, inclined air showers can be detected very 
well with the radio technique |17j . 

Some important questions have yet to be an- 
swered: Are the polarisation characteristics of 
the emission and the detailed angular correla- 
tions as expected by theory JHI? Results from 
the analysis of inclined air showers ^2] are con- 
sistent with the predicted signal polarisation, but 
only dedicated dual-polarisation measurements as 
are currently being prepared with LOPES will 
answer these questions definitely. What is the 
absolute strength of the emission, and is it con- 
sistent with the values predicted by theory |1U|? 
The LOPES30 calibrated measurements will an- 
swer this long-standing question soon. What is 
the scale parameter of the radio emission's ex- 
ponential radial decay, and how does it relate to 
the air shower geometry? LOPES30 will allow 
a much more detailed investigation of the per- 
event lateral radio profile than LOPES10. And 
finally, is the radio emission sensitive to the pri- 
mary mass of the EAS 20] and how can this be 
exploited? This question and many others can 
be investigated with the new CORSIKA-based 
Monte Carlo code [T2"] . 

4. The way to high energies 

The radio technique has high potential for ap- 
plication on large scales, and, consequently, at 



high energies. In particular, its ability to measure 
with nearly 100% duty cycle is a strong asset for 
combining it with surface detector arrays to take 
advantage of hybrid detection of EAS. One ap- 
plication of the radio technique on a much larger 
scale than LOPES or CODALEMA will be the 
use of LOFAR to study cosmic rays. 

Another vision is to use the radio technique for 
hybrid detection of ultra-high energy cosmic rays 
in the framework of the Pierre Auger Observatory 
(PAO). It is clear, however, that an application 
on this scale has requirements very different from 
those of already existing experiments. For once, 
it will be necessary to distribute the antennas on 
a grid much larger than that of the current exper- 
iments. A grid spacing of order 500 m or more is 
probably needed to make a large-scale radio ar- 
ray cost-effective. As a direct consequence, the 
organisation of such a radio array will have to be 
decentralized with mostly autonomous detectors 
(low power-consumption, ...) and wireless data 
transfer. The possibility to self-trigger on the ra- 
dio signal only would also be highly desirable (al- 
though not necessary in case of the PAO), and the 
feasibility of this is related to the grid spacing and 
the environmental noise levels (which luckily are 
much lower in Argentina than in environments 
such as the Forschungszentrum Karlsruhe) . 

A number of groups with expertise in radio 
detection of EAS are currently developing the 
necessary techniques, evaluating alternative an- 
tenna designs, and developing electronics and 
data links. Some small test cells are going to be 
set up inside the Pierre Auger Observatory start- 
ing September 2006 in order to evaluate their per- 
formance in the Argentinian noise environments. 

The experience and results gained — and to 
be gained — from LOPES at energies up to ap- 
prox. 10 18 eV are invaluable in their own right 
for the adaptation of the radio technique to large 
scales. In addition, LOPES has developed a new 
kind of self-triggering array called LOPES STAR . 
It was conceived to meet the requirements of a 
large scale decentralized array and is using a dif- 
ferent type of electronics and antennas. Detectors 
of the LOPES STAR type are going to be installed 
as part of the test cells in the Pierre Auger Obser- 
vatory. At the same time, they will be operated 



8 



and further developed in the environment of the 
original LOPES array, which continues to provide 
a unique test bed for the study of radio emission 
from EAS. 

5. Conclusions 

In the last few years, radio detection of cos- 
mic ray air showers has once again become a 
very active field of research. LOPES, as one of 
the projects studying radio emission from EAS 
with modern digital technology, has made impor- 
tant contributions. Radio emission from EAS is 
proven to be coherent at 40-80 MHz, and the 
field strength correlation with geomagnetic an- 
gle strongly indicates that the radiation is of 
dominantly geomagnetic origin. Aspects such as 
the absolute field strength, the detailed angu- 
lar correlations and the scale radius of the lat- 
eral dependence will soon be analysed with the 
new LOPES30 data. Our next-generation Monte 
Carlo code calculating the radio emission based 
on highly realistic CORSIKA-based air showers 
will also allow many new studies with unprece- 
dented detail. Unlike 40 years ago, the study of 
radio emission today is making good progress. 

A possible application of the radio technique 
on large scales for hybrid detection of ultra- 
high energy cosmic rays is now being studied 
in the framework of the Pierre Auger Observa- 
tory. These activities are still in a relatively 
early phase, but with the experience gained in 
the LOPES and CODALEMA experiments, the 
challenges to unlock the great potential of the ra- 
dio technique on large scales can be tackled. 
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